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Abstract. Cosmological weak lensing by the large scale structure of the Universe, cosmic shear, is coming of age 
as a powerful probe of the parameters describing the cosmological model and matter power spectrum. It comple- 
, ments Cosmic Microwave Background studies, by breaking degeneracies and providing a cross-check. Furthermore, 

^ ' upcoming cosmic shear surveys with photometric redshift information will enable the evolution of dark matter to 

fS) ' be studied, and even a crude separation of sources into redshift bins leads to improved constraints on parameters. 

CO . An important measure of the cosmic shear signal are the shear correlation functions; these can be directly calcu- 

' lated from data, and compared with theoretical expectations for different cosmological models and matter power 

spectra. We present a Monte Carlo method to quickly simulate mock cosmic shear surveys. One application of 
this method is in the determination of the full covariance matrix for the correlation functions; this includes red- 
shift binning and is applicable to arbitrary survey geometries. Terms arising from shot noise and cosmic variance 
(dominant on small and large scales respectively) are accounted for naturally. As an illustration of the use of such 
• covariance matrices, we consider to what degree confidence regions on parameters are tightened when redshift 

Oh' binning is employed. The parameters considered are those commonly discussed in cosmic shear analyses - the mat- 

ter density parameter f2m, dark energy density parameter (classical cosmological constant) JIa, power spectrum 
5^ ' normalisation as and shape parameter F. We incorporate our covariance matrices into a likelihood treatment, 

^ I and also use the Fisher formalism to explore a larger region of parameter space. Parameter uncertainties can be 

. decreased by a factor of ~ 4 — 8 (~ 5 — 10) with 2 (4) redshift bins. 
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1. Introduction obtained. For example, the CFHT Legacy survey 

. . ^ , ,. , . r ■,■ . • (http://www.cfht.Hawaii.edu/Science/CFHLS) will cover 

1 he statistics 01 the distorted images or distant galaxies, i-rr, j :^ • r- , > i u 4. u 

f , , , n 1 1 r deg"" m 5 optical bands, with a smaller area to be 

gravitationally lensed by the tidal gravitational held or , j ■ t j ta u j 

f . . ,,r observed m J and K bands, 
intervening matter inhomogeneities, contain a wealth of 

information about the power spectrum of the dark and o^^er to compare these observations with predic- 
luminous matter in the Universe, and the underlying cos- ^ions for various cosmological models and matter power 
mological parameters. The importance of "cosmic shear" spectra, different two-point statistics of galaxy ellipticities 
as a cosmological tool was proposed in the early 1990s by can be employed, all of which are filtered versions of the 
Blandford et al. (1991), Miralda-Escude (1991) and Kaiser convergence power spectrum. Here, we focus on the gravi- 
(1992). Further analytic and numerical work (e.g.. Kaiser tational shear correlation functions, which can be directly 
1998; Schneider et al. 1998; White & Hu 2000) took into obtained from the data as described in Sect.EJ 
account the increased power on small scales, resulting from This quest for the parameters describing the matter 
the non-linear evolution of the power spectrum (Hamilton content and geometry of the Universe is limited by sev- 
et al. 1991; Peacock & Dodds 1996). eral sources of error, dominated by the dispersion in the 
The feasibility of cosmic shear studies was demon- intrinsic ellipticities of galaxies and by cosmic (sampling) 
strated in 2000, when four teams announced the first variance. The covariance (error) matrix is thus an essen- 
observational detections (Bacon et al. 2000; Kaiser tial ingredient in the extraction of parameters from data, 
et al. 2000; van Waerbeke et al. 2000; Wittman et or in parameter error estimate predictions. Schneider et 
al. 2000). Upcoming surveys will cover much larger al. (2002a) provide analytical approximations for the con- 
areas, and multicolour observations will enable pho- tributions to the covariance matrix. They consider the 
tometric redshift estimates for the galaxies to be case when the mean redshift of the population is known, 
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and data taken in a single contiguous area. Kilbinger & 
Schneider (2004) use a numerical approach to investigate 
the impact of survey geometry on parameter constraint. 
Using a Fisher matrix approach, which provides a lower- 
bound estimate of covariance, Hu (1999) has shown that 
even crude redshift information enables much tighter con- 
straints to be placed on cosmological parameters, com- 
pared with the case when only the mean redshift of the 
population is known. This study concentrated on the con- 
vergence power spectrum as the vehicle of cosmological 
information. 

Motivated by these studies, in this paper we demon- 
strate how numerical simulations can be used to estimate 
the full covariance matrix for the shear correlation func- 
tions in the presence of redshift information, and for ar- 
bitrary survey geometries. We consider auto- and cross- 
correlations for redshift bins (as in Hu 1999) and in addi- 
tion allow for cross-correlations between measurements of 
the shear signal at different angular scales. With covari- 
ance matrices in hand, we then investigate the improve- 
ment in parameter estimates due to redshift binning. 

Further details and derivations of the equations rele- 
vant to cosmic shear and weak lensing can be found in 
Bartelmann & Schneider (2001). For a recent review of 
cosmic shear and future prospects see van Waerbeke & 
Mellier (2003). 



2. Power spectrum and correlation functions 

Access to cosmological parameters is provided through the 
observable two-point statistics of the ellipticities of distant 
galaxies. In this section, we describe how these are related 
to the matter power spectrum, and to the underlying den- 
sity field. 



2.1. The convergence power spectrum 

The power spectrum Pk,{£) of the effective convergence, 
or equivalently of the shear P-y{£) (e.g. Bartelmann & 
Schneider 2001), is related to that of the density fluctu- 
ations, Ps{i), through a variant of Limber's equation in 
Fourier space (Kaiser 1998) 
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(1) 
(2) 
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where £ is the angular wave- vector, Fourier space con- 
jugate to 6. w is the comoving radial distance, K the 
curvature parameter. A value K = corresponds to 
rim + = 1, where Qa is the cosmological constant and 
J7ni the matter density parameter. The function W{w) ac- 
counts for the sources being distributed in redshift, where 



p{w')dw' is the comoving distance probability distribu- 
tion for the sources, a (w) is the scale factor normalised to 
a (ui = 0) = 1 and Ho is the Hubble constant. 

Splitting up the weak lensing survey in redshift, as 
in FigQ] defines a set of effective convergence and shear 
maps instead of a single one, giving more information on 
the evolution of the dark matter fluctuations and therefore 
enabling tighter constraints to be placed on cosmological 
parameters. Auto- and cross-correlation of these maps in- 
troduce a whole set of power spectra, generalising Eq. : 

^ ^ ' " 4c4 
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with p^'^^ {w) being the normalised distribution in comov- 
ing distance inside the ith bin, where i runs between 
1 and the number of redshift bins N^. Pi"'* are auto- 
correlation power spectra, whereas Pi*"''' with i ^ j are 
cross-correlation power spectra. 



2.2. Shear correlation functions 

Constraints can be placed on cosmological parameters 
using the directly observable shear correlation functions, 
which we now turn to. 

The basis which underpins the use of the distorted im- 
ages of distant galaxies in weak lensing studies is a trans- 
formation relating the source, e*-''-', and image, e, (complex) 
ellipticities to the tidal gravitational field of density inho- 
mogeneities (for definitions see Bartelmann & Schneider 
2001). We focus on the non-critical regime where 



1 + .g*e(") 



(5) 



where g = 7/(1 — k) is the reduced shear. 

Empirically the probability distribution function (pdf ) 
of the galaxies' intrinsic ellipticities is a truncated 
Gaussian for both the real and imaginary parts of e^*'-': 



exp (~|6(-)|Va^,.,) 
Tia%, [l~exp(-l/(T2^,)] ' 



(6) 



where (Tj(a) is the intrinsic ellipticity dispersion. 

As in Schneider et al. (2002b), the shear correlation 
functions are defined as 



C±(^) = (7t7t) ± (7x7x> 



'2^ 



Jo,4(^^) P.{t) , (7) 



where J„ are n-ih order Bessel functions of the first kind; 
7t and 7x are the tangential and cross shear components 
respectively. From now on we focus on ^_|_, since this con- 
tains most of the cosmological information on the scales 
of interest. 
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2.3. Choice of cosmology and matter power spectrum 

Unless otherwise stated, our cosmology throughout is a 
ACDM model with fl^ = 0.3, f^A = 0.7 and Hq = 
70kms^^ Mpc~^. A scale-invariant {n — 1, Harrison- 
Zel'dovich) spectrum of primordial fluctuations is as- 
sumed. Predicting the shear correlation functions requires 
a model for the redshift evolution of the 3-D power spec- 
trum. We use the fitting formula of Bardeen et al. (1986; 
BBKS) for the transfer function, and the Peacock and 
Dodds (1996) prescription for evolution in the nonlin- 
ear regime. The power spectrum normalisation is parame- 
terised with as = 0.9, and T = 0.21. Quantities calculated 
for this fiducial cosmology/power spectrum will be super- 
scripted with a "t" . 

3. Simulating cosmic shear surveys 

In this section we describe the method we used to make 
Monte-Carlo simulations of cosmic shear surveys. The im- 
plementation makes the simulations as computationally 
inexpensive as possible, i.e. without invoking N-body sim- 
ulations. 

Calculating the Icnsing signal by ray tracing through 
N-body simulations has become a common tool for making 
simulated weak lensing surveys (see e.g. Blandford et al. 
1991; Wambsganss, Cen & Ostriker 1998; Jain et al. 2000). 

We take a different path here, because only the two- 
point statistics of weak lensing is considered. This allows 
us to reduce the computational effort by expressing the 
fields of the shear and convergence as random Gaussian 
fields having the same power spectrum as the correspond- 
ing fields from the N-body approach. In the weak lens- 
ing regime, ray-tracing is well described by the Born ap- 
proximation which ignores the effects of lens-lens coupling 
and deviations of light rays from the fiducial path (see 
White & Hu 2000). The task of calculating the required 
power spectra then becomes relatively straightforward, be- 
cause these can be shown to be linear functions of the 
three-dimensional evolving dark matter power spectrum. 
The accuracy of the results depends on how accurately 
that three-dimensional power spectrum is known. 

Since we consider the two-point cosmic shear statistics 
in an area of relatively small angular size, we can represent 
the cosmic shear fields by random Gaussian fields in a flat 
sky approximation. 

The practical advantage of simulating a single random 
Gaussian field <5 (r) - a homogeneous and isotropic random 
Gaussian to be exact - is the fact that from the Fourier 
coefficients 

Ck^^ J^drS (r) exp (ik • r) (8) 

of such a random field only a pair of coefficients is corre- 
lated 

(ckC-k') = -^(^kk'-Pk , (9) 

where (5d is the Dirac delta function and Pk is the power 
spectrum of the random field. The volume V of the 
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Gaussian field is in our case simply the area on the sky 
covered by the field. 

In our work, however, the situation is a bit more com- 
plicated than that: the ellipticities of galaxies belonging to 
different redshift bins are correlated as well as the elliptic- 
ities of galaxies at different angular positions. Therefore, 
when more than one redshift bin is considered, several 
Gaussian fields - cosmic shear maps - with prescribed 
cross-correlations have to be simulated simultaneously. A 
way to do this in general on a regular grid for real Gaussian 
fields is shown in the following subsection. The subsection 
thereafter explains how we used this approach for simu- 
lating mock cosmic shear surveys. 

3.1. Realisations of Correlated Gaussian Fields 

According to condition the pair Ck and c_k is corre- 
lated. In this section we restrict ourselves to real Gaussian 
fields with 6 (r) = S* (r) . This introduces an additional 
condition that follows from the definition ||SJ) of the Ck: 

Ck = c*_k . (10) 
In particular, for real Gaussian fields we thereby have 

(CkCk'> = -^-Pk-^kk' , (11) 

<Jkk' being the Kronecker symbol. The conditions © and 
(|10|l are easily accounted for if, say, only the Ck for half 
of the spatial frequencies k are worked out and the c_k 
frequencies are set accordingly. Hence, for our choice, if we 
talk about Ck we actually mean only Fourier coefficients 
in the right half-plane. 

Furthermore, the real and imaginary parts of Ck are 
uncorrelated, and both follow the same Gaussian pdf. This 
pdf has zero mean^ and a variance dk that is expressed in 
terms of the power spectrum Pk describing the two-point 
correlations of the fluctuations in the Gaussian field (see 
e.g. Peacock 2001) 

= ^Pk . (12) 

The procedure for making one Gaussian field realisa- 
tion requires two steps: 1. drawing numbers for the real 
and imaginary parts for every Ck with a Gaussian random 
number generator, and 2. transformation of this Fourier 
space representation to real space in order to obtain the 
field realisation. For the second step we used an FFT al- 
gorithm from Press et al. (1992)^. 

This procedure also holds when realisations of more 
than one, but uncorrelated Gaussian fields are desired. 

^ In the case that S (r) has a non zero mean, (ck) for k = 
becomes different from zero. 

^ As in FFT the matrix of the Fourier coefficients contains 
Ck that share the same matrix elements with c_k, one has for 
these particular coefficients to set the imaginary parts to zero 
and to increase ak by the factor \/2. The latter is necessary to 
guarantee that the variance of tire modulus of Ck is still correct. 
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"Uncorrelated" means that if we denote the Fourier co- 
efficients of, say, N Gaussian random fields by cj^-' with 
i = 1..N then we expect for those fields the relation 



1 



(13) 



where 6ij is also a Kronecker symbol, this time for the 
Gaussian field indices. P^"^ is the previously introduced 
power spectrum, or auto- correlation power spectrum, of 
the ith random field. Thus, here correlations between Ck 
of different random fields vanish. 

For the purposes of this work, however, we need to 
be able to allow for cross-correlations between different 
random fields i 7^ j in a defined manner, like 



P^''^ is for i j the cross- correlation power spectrum. 
Like for the auto-correlations, only certain pairs of Fourier 
coefficients of different Gaussian fields are correlated. This 
follows from the assumption that the cross-correlations are 



.0) 



1 



(14) 



(ij) 



P, 



homogeneous, too. Note that Pj 

In order to find a recipe for making realisations of that 
kind, we make the Ansatz that the N Fourier coefficients 
c"^ are a linear transformation Ak of N different uncor- 
related coefficients df^ with an equal Gaussian pdf for the 
real and imaginary parts, zero mean and a l/\/2 disper- 



(9) 



(15) 



Ak is a iV X iV linear transformation matrix. The linearity 
of this transformation accounts for the fact that the result- 
ing set of coefficients cj^*'' still obeys a Gaussian statistics, 
because linear combinations of Gaussian random variables 
are also Gaussian. 

Since for real Gaussian fields real and imaginary parts 
of the Fourier coefficients are not correlated, and by our 
Ansatz neither are the real and imaginary parts of d!^\ 
only real numbers for the components [Ak]*^ are allowed; 
an additional imaginary part of Ak would mix real and 
imaginary parts of dj^*-* , thereby possibly introducing cor- 
relations between real and imaginary parts in Ck. A ma- 
trix Ak that is purely imaginary would be an alternative 
choice, though. 

Eqs. 115(1 can now be combined to give 



q,r 

5][Ak]\[Akr,v 



Ar) 



T]9 
k -i 



(16) 



where A^ denotes the transpose of Ak. Hence, together 
with equation (|14|1 this puts another constraint on the 
matrix Ak, namely 



1 piv) 



(17) 



For convenience we introduce the power matrix defined as 
[PkYj = yP^''^ to abbreviate this equation: 



AkAj 



(18) 

The power matrix is the covariance matrix between the 
Fourier coefficients of a set of Gaussian fields for a certain 
k. 

This shorthand of N"^ equations does not uniquely de- 
termine the matrix Ak, because it contains only N{N -\- 
l)/2 linearly independent equations, since both the matrix 
on the Ihs and the matrix product on the rhs are symmet- 
ric. As there are no further constraints on Ak, we are al- 
lowed to set the remaining N'^-N{N-\-l)/2 ^ N{N-l)/2 
constraints of Ak as we like. We do this by assuming that 
Ak is symmetric, so that we finally obtain 

Pk = A2 ^ Ak = Vl\. (19) 

In general the square root is not unique (see e.g. 
Higham 1997). However, we are already satisfied with one 
particular solution to this problem. In order to determine 
such a solution, note that Pk is a symmetric positive 
(semi)definite matrix, which is ensured by the properties 
of the power spectra the power matrix consists of: 



< Pi 'P, 



(") 



(20) 
(21) 



Therefore, this matrix can uniquely be decomposed into 

Pk = Rl Dk Rk , (22) 

where Rk is an orthogonal matrix whose column vectors 
are the eigenvectors of Pk, while their corresponding, al- 
ways real and positive, eigenvalues Xi are on the diagonal 
of the diagonal matrix Dk = diag (Ai, A2, An). As one 
particular square root we pick out 



(23) 




which is a solution due to 



/D];^/DJ; = Dk. 

To sum up, for every k mode considered, the process 
for the realisation of correlated Gaussian random fields 
requires one to find the square root Ak of the power ma- 
trix Pk. This defines a linear transformation for a vec- 
tor of uncorrelated random complex numbers (real and 
imaginary part of the same coefficient are uncorrelated, 
too) with zero mean, real and imaginary parts obeying a 
Gaussian pdf with 1/V2 variance. Applying Ak yields a 
vector of Fourier coefficients belonging to the realisations 
of the correlated Gaussian random fields. Due to pU|) this 
is performed only for one half of the spatial frequencies 
considered. The other half is set accordingly to fulfil this 
condition. 
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3.2. Simulating the weak leasing survey 




Fig. 1. Galaxies are binned together according to their 
redshift, the boundaries of the pairwise adjacent redshift 
bins are Wi with i = O...Nz (here as an example Nz — 4). 
For every redshift bin the reduced shear field is calculated, 
averaging over the redshift distribution inside the bin. 



Each galaxy in the mock galaxy catalogue is defined by 
an angular position, an ellipticity e and a redshift bin it 
belongs to. The ellipticity of the isophotes of a galaxy is 
determined by the intrinsic shape of a galaxy ef"^ and the 
reduced shear g at the position of the galaxy (see Eq.EJ. 
The reduced shear is a function of the convergence k and 
shear 7 which have to be simulated for each redshift bin as 
a map covering the simulated area. Here we assume that 
the galaxies are binned into Nz pairwise adjacent redshift 
bins, chopping off the redshift distribution 



p{z) 



1 



1 



z3r(3//3)/3 



exp 



(24) 



as in Fig.^ This empirical distribution with /3 = 1.5 and 
zo = 1-0 is based on deep field surveys (see e.g. Small et 
al. 1995). The total number of galaxies inside the field, 
with chosen size of 5° x 5°, is set to be « 2.7 x 10^, to 
get an average of 30 galaxies per arcmin^. Moreover, the 
galaxies are assumed to be randomly distributed over the 
field of view. 

The method of the last subsection is used to work 
out the convergence maps in Fourier space on a grid of 
2048 X 2048 pixels for all redshift bins. As input the power 
matrix Pk, consisting of the auto- and cross-correlation 
power spectra of these convergence maps, specified by the 
Eqs. (gj, is needed. 

In the next step, the shear maps are obtained from the 
convergence maps using the relation 



7-6 



(■2 ~ 



t = 



7£ and are the Fourier coefficients of the shear and con- 
vergence fields for the angular frequency £, respectively. 
This relation stems from the fact that both shear and 
convergence are linearly related to a potential function. 
For every galaxy, shear and convergence are then com- 
bined with the intrinsic ellipticity , randomly drawn 
from the pdf Eq. © using a^f^ — 0.3, to compute the final 
ellipticity of the galaxy via Eq. 

Both angular size A and number of pixels Np along 
one axis - the sampling size - limits the number of fluc- 
tuation modes accounted for in the simulated data. This 
means, since we are lacking fluctuations on scales outside 
of b°/ {2Np) < e < 5°, equivalent to < (. < Cax, 
that we have less correlation in the cosmic shear fields 
than expected (Eq.[7Jl 



2,71 



V2A ' 



(26) 
(27) 



(25) 



The values for the limits are estimates for a square field; 
the limits are not clearly defined, because the number of 
^-modes in the EFT matrix becomes very small near the 
cutoffs. One solution to this problem is to artificially set 
a clearly defined range within the interval [^mim^max], or, 
as we have done, to find a best fit cutoff. This is found by 
varying the cutoffs to obtain closest agreement between 
the theoretical two-point correlation and the ensemble av- 
erage of all Monte-Carlo realisations. 

In total, we simulated two data sets. The first data set 
consists of Nf = 795 independent realisations each 5° x 5°. 
The redshift distribution of the galaxies was split into 2 
bins at a redshift cut Zcut = 1-25, and the distribution is 
truncated at z = 3. For our fiducial surveys, we randomly 
selected 10 sub-fields, each of 1.25° x 1.25°, from differ- 
ent large realisations. This was done for two reasons: 1. to 
reduce the computation time, since for 10 shear maps we 
require only one realisation, and 2. sub-fields are less af- 
fected by ^ijiin that necessarily enters the simulations due 
to the finite realisation area. 

The second data set has 4 redshift bins, with Zcut = 
0.75,1.5,2.25,3.0 where the last value is the truncation 
redshift. It has Nf — 266 independent realisations. The 
fiducial surveys from this data set consist of single sub- 
fields of size 1.25° x 1.25°. 

For both data sets, ^+ was estimated (see next section) 
for Nao = 65 angular separation bins, ranging from about 
2'.0 to 40'.0. For the first (second) data set, the correlation 
functions were subsequently averaged for 10 (1) sub-fields 
in order to simulate cosmic shear surveys consisting of 10 
(1) independent data fields. 

In a further step, the cross- and auto-correlation of 
the cosmic shear between the shear maps were, according 
to appendix ^ combined to yield the cosmic shear cor- 
relations for a coarser redshift binning; in each step the 
number of redshift bins was reduced by one by combining 
two neighbouring bins. This process gave for the first data 
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Table 1. The final column denotes the name given to a 
particular binning of data. The entries in the columns Zi 
show the corresponding cuts in redshift, Zcut- 





Z\ 


Z2 


Z3 


Z4 


Name 





0.75 


1.5 


2.25 


3.0 


4bins 





1.5 


2.25 


3.0 




3binsl 





0.75 


2.25 


3.0 




3binsII 





0.75 


1.5 


3.0 




3binsIII 





2.25 


3.0 






2binsl 





1.5 


3.0 






2binsII 





0.75 


3.0 






2binsIII 





1.25 


3.0 






2binsIV 





3.0 








onebin 



set, apart from the original data, the shear correlation of 
one redshift bin with boundaries 2 = and z = 3. The 
second data set allows more freedom of choice for combin- 
ing redshift bins, so that we are able to construct several 
data sets with three and two redshift bins. TableQllists the 
different redshift bins and reference names, all extracted 
from the two original data sets. 

Below, this data is used to study the improvement in 
the statistical uncertainties of the cosmological parameter 
estimates, if one has more information on the redshifts of 
the galaxies. 

4. Estimating ^+ 

To estimate the two-point correlator ^+ between the 
galaxy ellipticities - depending on position 0^ and red- 
shift bin - inside of the sub-fields we use the estimator 



1 



X ^ WiWj {eitejt + ejxejx) Ae {\@^ - 

ij 



_/ 1 for 61 - A0/2 <(l)<e + Ae/2 







otherwise 



(28) 



as mentioned in Schneider et al. (2002a), with A^ being 
the width of the angular bins. Since we are dealing with 
simulated data here, there is no need to weight galaxies 
with respect to their ellipticity. Therefore, we set Wi = 1 
for every galaxy. 

Although mathematically simple, it takes quite a time 
to evaluate the estimator due to the large number of 
galaxy pairs. To speed up the whole procedure we put 
a grid of rectangular cells of size A9 x A6 over the sub- 
field in question and compute the number Nij of galaxies 
and the mean of their ellipticities e^- inside every cell. The 
index ij indicates the position of the cell inside the grid. 
This means we are representing galaxies inside the same 
cell by a single data point with weight Nij and ellipticity 



lij. In particular, all galaxies inside this cell are assumed 
to be placed at the same position. The estimator of ^+ for 
this rearranged data set can be shown to be 



^ijx^klx ) Ag 



&kl\) 



ij,kl 



Np {9) = N^j^kiAe (|0„- - @ki\) , (29) 

ij.kl 

where 0^^ represents the angular position of cell ij. 

The advantage of this approach is obvious: instead of 
considering iV^ pairs {N is the number of galaxies) we 
have to consider only N"^ pairs, where Nc is the number 
of grid cells. Thus, the number of pairs depends only on 
the cell size and not on the number of galaxies. Hence, this 
method pays off once the cell size becomes large enough, 
making the number of cells smaller than the number of 
galaxies. Moreover, in order to find all galaxies at some 
distance from a certain cell we no longer have to check 
all galaxies, but only neighbouring cells which are easy to 
find by the grid index. 

The approach becomes inaccurate, however, for small 
angular bins, because for these the assumption that cell- 
galaxies are essentially concentrated into one single point 
is particularly inaccurate. By comparing the ensemble av- 
erage of ^+ with the theoretical ^+ we find that after the 
third angular bin this approximation becomes accurate 
enough. For our purposes, this approach is completely suf- 
ficient. A better and more sophisticated approach can be 
found in Pen & Zhang (2003). 

For the case with Nz = 2, with the division at Zcut = 
1.25 (2binsIV) Fig.01 shows the close agreement between 
the correlation and cross-correlation functions, averaged 
over 7950 sub-fields, with the analytical prediction for the 
fiducial ACDM cosmological model, obtained from Eq. 0. 
Shown are comparisons for the lower (L) and upper (U) 
redshift bins, and cross-correlation (LU). To account for 
finite field size in our numerical work, Imi-a = 27r/14.9° 
in the integration. As noted above, since there is no well- 
defined cut-off, this value of ^min was determined by allow- 
ing it to vary while performing a least-squares fit of y uj 

to ^Cl,u,lu^: so obtaining the inverse variance weighted 
mean fmin- A cut-off at high ^ is not critical since in this 
regime the power-spectrum amplitude is much lower. 

5. Estimating the covariance of ^± 

We now outline how the covariance matrix of ^± is es- 
timated, for the case of Nz = 2 redshift bins, with the 
division at z^nt = 1-25 (2binsIV). As described above, our 
mock survey consists of 10 uncorrelated fields. An angle 
bracket denotes averaging over all 7950 sub-fields. Note 
that we may drop the L, U and LU sub-scripts for ease of 
notation. 

When no redshift binning is assumed, it is computa- 
tionally advantageous to determine the shear correlation 
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LU 



U 



LU 



U 



Fig. 2. Schematic illustration of the symmetric covari- 
ance matrix C for the case where there are Nz = 2 
source redshift bins and Nab = 5 angular separation bins. 
Combinations of ^l.u,lu identify covariance terms of the 
form given in Eq. (|31|l . 

function by combining those determined for the case with 
redshift binning: 



C = "l^l + 2nL7iu'^LU 



(30) 



where rtL,u are the fraction of sources in the lower and 
upper bins respectively. A more general relation between 
^ and the cross- and auto-correlations of the shear from 
more than two redshift bins can be found in appendix [XI 
The impact of cosmic variance is taken into account by 
determining correlation functions for each of the Nf — 795 
independent surveys. The covariance matrix between bin 
i and j is determined using 



(31) 



N, 



where the outer average is performed over the Nf — 
795 surveys. If redshift binning is considered, there are 
-^z {^z + 1) /2 combinations of correlation and cross- 
correlation functions and hence C is a matrix composed 
of [Nz {Nz + 1) /2f blocks. Fig.H illustrates this for the 
simplified case where Nz = 2 and Nab = 5; for exam- 
ple the block in the upper left of the matrix corresponds 

to elements = - (c)^^) - ™^ 

the shaded entry to C2,3. The block in the middle row, 
left column, corresponds to covariance elements between 
the cross-correlation and the lower redshift bin, with the 
shaded entry being C8,2- The bins denoted by i and j ex- 
tend over A0 bins, repeated for each redshift auto- and 
cross-correlation bin. 



^11 ^12 ^13 ^14 ^2 ^3 ^4 ^33 ^34 ^44 




Fig. 3. The covariance matrix C determined from our sim- 
ulations, for Nz = 4, Nab = 65 and a survey consisting 
of 1 sub-field. Different blocks correspond to auto- and 
cross-correlations between the redshift bins. Inside these 
blocks are auto-and cross-correlations for angular separa- 
tion bins. 
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T(5 



10 

6 (arcmin) 

Fig. 4. Comparison of the analytical (lines) and numeri- 
cal (symbols) shear correlation and cross-correlation func- 
tions ^L,LU,u (lower, middle and upper lines/symbols). 



A representation of the covariance matrix determined 
from our simulations with TV^ = 4 is shown in Fig.|3| Note 
that C has a strong diagonal, although it is not strictly 
diagonally dominant. 

Our covariance matrix for the case of no redshift bin- 
ning is consistent (< 10% difference) with the treatment of 
Schneider et al. (2002a), and with Kilbinger & Schneider 
(2004) who adopted the same assumption of Gaussianity. 
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6. An application: Constraints on cosmological 
parameters 

The rather featureless two-point shear correlation function 
^,+ (0) or corresponding convergence (shear) power spec- 
trum P^ii) leads to strong degeneracies amongst the pa- 
rameters that can be derived from cosmic shear surveys. 
An indication of the degree of degeneracy is the behaviour 
of the partial derivatives of ^+ with respect to each param- 
eter TTi (see King & Schneider 2003 for such a comparison) , 
or using a Fisher matrix analysis as in Sect. 16. 'ill 

External sources of information often provide compli- 
mentary constraints: for example, confidence regions in 
the Jim — as plane derived from weak lensing are almost 
orthogonal to those from the analysis of CMB data (e.g. 
van Waerbeke et al. 2002), lifting this well known degen- 
eracy (e.g. Bernardeau, van Waerbeke & Mellier 1997). 

In this section we consider the extent to which crude 
redshift information for sources used in a lensing analy- 
sis decreases the expected errors in the fl^ — erg , — T 
and (78 — r planes. Since we are interested in the influ- 
ence of redshift binning on parameter degeneracies, hid- 
den parameters are assumed to be perfectly known. As 
described above, we focus on the information provided 
by the shear two-point correlation function ^_|- (we may 
drop the subscript). We restrict this application to 
the case oi = 2 (2binsIV). A larger parameter space 
is then explored using the covariance matrix derived from 
simulations in a Fisher analysis for the cases — 2,3,4. 



6.1. Obtaining confidence regions in the 
— r and as — T planes 



0-8, 



We now determine and compare the likelihood contours in 
the f2„i — (Tg , ~ r and cg — F planes for the cases with 
and without redshift binning. The likelihood function is 
given by 



1 



(27r)«/2|C|i/2 
1 



(32) 



]^exp 



where n is the number of rows (or columns) of the covari- 
ance matrix C and ^(tt) are theoretical correlation func- 
tions determined on a grid in parameter space. 

The log-likelihood function is distributed as so 
that 



X 



(33) 



Confidence contours can be drawn in this x^-surface, rel- 
ative to the minimum (zero) at f (tt) = In Figs. I5I7I the 
confidence contours are shown for each of the Qj^i — cg, 
— F and F — dg planes, with and without redshift bin- 
ning. Note that while is varied, we keep + = 1. 
To highlight the difference and avoid confusion, we plot 
contours for a single value of Ax^. 



1.3 
1.1 
Og 0.9 
0.7 
0.5 



0.1 0.2 0.3 



0.4 



0.5 0.6 0.7 



Fig. 5. Expected constraints in the flm ~ erg plane plotted 
for Ax^ = 4.61 (90% confidence) with (inner contour) 
and without (outer contour) redshift binning. The survey 
consists of 10 square uncorrelated sub-fields, each 1.25° on 
a side. The redshift distribution and binning are described 
in the text. 



m 




Fig. 6. Expected constraints in the F — ilm plane plotted 
for Ax^ = 4.61 (90% confidence) with (inner contour) 
and without (outer contour) redshift binning. The survey 
is the same as in Fig.|21 



6.2. Fisher information 

The Fisher matrix (Fisher 1935) gives a handle on the 
question as to how accurately model parameters can be 
estimated from a given data set. In this section, we will 
use this method to examine quantitatively the increase of 
information on the cosmological parameters f^m, JIa, erg 
and F when the number of redshift bins, thus the knowl- 
edge of the three-dimensional distribution of the galaxies, 
is increased. Note that we no longer impose the condition 
f2iji -|- Q\ = 1. After a brief introduction to this topic we 
apply the Fisher statistics to our simulated data. 

6.2.1. Fisher Formalism 

In general, one uses data points from a measurement to 
infer model parameters tt^ based on a theoretical model. 



Simon, King & Schneider: 3D correlations of cosmic shear 



9 



Oo 0.9 




Fig. 7. Expected constraints in the T — as plane plotted 
for Ax^ = 4.61 (90% confidence) with (inner contour) 
and without (outer contour) redshift binning. The survey 
is the same as in Fig.O 

As the measurements are polluted by noise, we cannot 
expect to exactly obtain the data points ^ (tt) predicted by 
our model. But we can try to find a combination of model 
parameters TTi that predict data points as close as possible 
to the actual measurement. The closeness is decided on 
the grounds of a statistical estimator. The covariance of 
the parameter uncertainties 



(34) 



1/2 



with A-Ki = (^{tt^) — (^i)^j is related to the so-called 
Fisher information matrix through 



dlTidTTj 



(35) 



£ corresponds to the likelihood for obtaining the mea- 
surement ^ keeping the underlying model parameter tt^ 
fixed. See for example Tegmark et al. (1997) and refer- 
ences therein for a more detailed description. 

It follows from statistics that the la scatter of the 
estimated parameters is (Cramer-Rao inequality) 



AtT,; > 



(36) 



where commonly the lower limit is taken to be the esti- 
mate for AiTi . To quantify the degeneracies in the param- 
eter estimate, we evaluate the correlation of the estimate's 
uncertainty contained in F: 



(37) 



as, for example, in Huterer (2002). Highly correlated 
or anti-correlated Att^ and Attj are called degenerate, 
whereas no correlation means no degeneracy (for the fidu- 
cial model) . Another piece of information that can be ex- 
tracted from the Fisher matrix is the orientation of the 
error ellipsoid in parameter space, which is defined by the 



eigenvectors of F. This corresponds to the directions of 
degeneracies. 

In the application of this formalism in Sect . 16. 2"^ be- 
low, we will look at situations where some of the model 
parameters are assumed to be known a priori. In this case, 
they are no longer free parameters that have to be esti- 
mated from measured data points, so that the size of the 
Fisher matrix reduces according to the number of param- 
eters fixed. This amounts to removing rows and columns 
from the general Fisher matrix, one for each fixed param- 
eter, so that these cases can be considered by simply look- 
ing at sub-matrices of the largest Fisher matrix. Taking 
all conceivable sub-matrices enables the exploration of all 
possible combinations of fixed (strong prior) and free pa- 
rameters. 

In practice, one uses the approximation (|32|l for the 
likelihood function £ (where i^* = ^), so that the Fisher 
information matrix is approximately 



kl 



diTi 



(38) 



which is exact in the case of pure Gaussian statistics, but 
may be used as a good approximation for the valley in 
parameter space in which the minimum of — log £ lies. 
Again, C is the covariance of the measured data points 
and ^ (tt) the vector of modelled data points in absence of 
noise. 

6.2.2. Application of the Fisher formalism 

Now we use the Fisher formalism to estimate constraints 
on various combinations of parameters, with different de- 
grees of redshift binning. First, we evaluate Eq. I|38|l using 
the covariance matrix from our fiducial survey consisting 
of 10 independent sub-fields, Nz = 2 (with Zcut — 1-25), 
and N^e = 65. The procedure is repeated for the covari- 
ance matrix for the coarser = 1 binning. Table|5]shows 
the percentage error for A^^ = 2 as opposed to Nz = 1 for 
the same set of free and fixed parameters. 

We extend the treatment to a larger number of redshift 
bins, in the context of the survey consisting of 1 sub-field. 
Again Eq. (|38|l is calculated, this time using the covariance 
matrices from the simulations with Nz = 4, and those 
from coarser binning {Nz — 3, 2, 1) of this data set. TableOl 
lists the errors for Nz — 4, 3, 2 as a percentage of the 
Nz — 1 error. 

In order to investigate the degeneracies of the param- 
eter estimates, we concentrate on the case that no priors 
are given. For this particular situation, the gain by intro- 
ducing redshift binning is largest (see Tables|2] and O . In 
Fig.lHlwe plot the correlations of the errors in the parame- 
ter estimates for different pairs of parameters and different 
numbers of redshift bins. If more than one redshift binning 
for the same number of redshift bins is available in our 
data set, we indicate the scatter of correlation coefficients 
by error bars. Some scatter indicates that the correlations 
can be changed slightly by varying the bin limits. The 
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strong correlation between the estimates of i7ni and is 
only marginally affected by redshift binning. This is also 
the case for fixed F and/or J^a (not shown). 



0.5 



-0.5 



-1 n 



■■X- 



12 3 4 

redshift bins 

Fig. 8. Correlations of the errors in the parameter esti- 
mates for different pairs of parameters and numbers of 
redshift bins as derived from the Fisher matrix; only the 
case with no fixed priors is considered. Error bars denote 
the variance in the correlations for the different redshift 
binnings for the same number of bins (only for 2 and 3 
bins). The data points are slightly shifted to avoid over- 
lapping. 



Table 2. Uncertainties in the parameter estimates ac- 
cording to the Fisher formalism, for our fiducial survey 
of 10 uncorrelated sub-fields. The first data set (2binsIV) 
with Nz — 2 and Zcut — 1-25 is used. Columns with dots 
"." denote fixed parameters (strong priors). Uncertainties 
in the top panel are absolute values for a single redshift 
bin. Those in the lower panel are for Nz = 2, quoted 
as a percentage of the single redshift bin {Nz = 1) case. 
For instance, with no fixed parameters, Afl\ = 0.26 with 
Nz = 2 (i.e. 13% of the Nz = 1 error). |F-i| denotes the 
determinant of the inverse of the Fisher matrix; its square 
root is proportional to the volume of the error ellipsoid in 
parameter space. The nth root, with n being the number 
of free parameters, defines a typical size of the error el- 
lipsoid; this size is proportional to the geometric mean of 
the lengths of the principal ellipsoid axes. 
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7. Discussion 

The average shear correlation functions obtained from 
our numerical simulations are in good agreement with 
those obtained analytically, as was illustrated in Fig.0] 
We also pointed out that their covariance is compatible 
with Schneider et al. (2002a) and Kilbinger & Schneider 
(2003). 

Our treatment is only strictly valid for Gaussian den- 
sity fields and is a good approximation for scales greater 
than ^ 10', giving a lower limit on the covariance at 
smaller scales (e.g. van Waerbeke et al. 2002). A more 
accurate covariance matrix is possible, though. According 
to Schneider et al. (2002a) (section 4 therein), the covari- 
ance matrix of ^_|_ may be decomposed into three terms 

= a^X.j + ajY.j + , (39) 

where X, Y and Z are some functions. X and Y are func- 
tions of the two-point correlation of cosmic shear and con- 
sequently insensitive to non-Gaussian features of the field. 
Z, however, depends linearly on the four-point correlation 
of cosmic shear which in Schneider et al. (2002a) is worked 
out by assuming a Gaussian field; this factorises Z into a 
sum of products of two-point correlators only. In the hi- 
erarchical clustering regime, the four-point correlation of 



0.02 


0.10 


0.02 


0.06 


0.15 
0.10 
0.02 
0.06 


13% 


13% 


17% 


26% 


47% 




52% 


57% 


64% 


72% 


53% 




52% 


49% 


73% 


43% 


38% 




68% 


65% 


32% 


24% 


46% 




57% 


50% 


45% 






64% 


86% 




85% 




88% 


89% 






81% 


85% 




65% 


71% 




76% 




80% 




81% 


81% 






87% 


80% 


84% 


85% 








85% 




81% 






81% 






89% 




89% 








82% 


82% 



the random field differs from that value only by a constant 
scale- independent factor Q, the so-called hierarchical am- 
plitude (see e.g. Bernardeau et al. 2002). Thus, hierarchi- 
cal clustering increases the component Z simply by the 
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Table 3. The notation is identical to Table |21 The second 
data set was used to obtain these values, with various bin- 
nings in redshift denoted by the final column, and with co- 
variance matrices calculated for 1 sub-field. Uncertainties 
are again quoted as a percentage of that for Nz = I. 
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factor Q. This could be included in dj by the following 
two steps: i) calculating dj = c|j^^ by setting the intrinsic 

(2) 

noise — 0, ii) recalculating = Cl^ this time with 
the intrinsic noise turned on. The final covariance matrix 
Cij, accounting for Q, is obtained by 

a,=C«(Q-l) + clf . (40) 

As an illustration of the use of numerically derived 
covariance matrices, we have considered to what degree 
redshift information tightens the confidence regions on 
cosmological parameters. Note that we required rather 
crude redshift information: photometric redshift estimates 
for sources can be obtained from multi-colour observa- 
tions, with typical accuracies in cr^ of 0.1 or better (e.g. 
Bolzonella et al. 2000). 



We have investigated the information contained in ^_|_, 
rather than the shear (or equivalently convergence) power 
spectrum P^, since it is directly obtained from the statis- 
tics of the distorted images of distant galaxies. Various 
estimators of have been proposed, requiring the spatial 
distribution of the shear (e.g. Kaiser 1998; Hu & White 
2001) or the shear correlation functions (Schneider et al. 
2002a). Note that a calculation similar to ours but using 
Pk would formally require one to use the full covariance 
matrix for the associated estimator. However, as noted in 
Schneider et al. (2002a), the band-power estimates for P^ 
do decorrelate rather quickly. 

The first data set consists of 2 redshift bins, and co- 
variance matrices estimated with and without binning for 
a survey with 10 uncorrelated (1.25° on a side) sub-fields 
(i.e. selected from different realisations). Constraints on 
pairs of cosmological parameters using a likelihood treat- 
ment, with and without redshift information, are shown 
in Figs . 15171 Since our goal here is to study the benefit of 
redshift information in cosmic shear studies, we do not 
adopt priors from WMAP or other probes of large scale 
structure which might confuse the issue. In both cases the 
redshift distribution is assumed to be known. Assigning 
sources to 2 redshift bins tightens the confidence regions 
in all cases. 

Noteworthy are the tightened upper limits on T in the 
r — (Tg and r^in — r planes when binning is employed. 
The constraints on F in both planes are rather asym- 
metric, with the confidence regions being more extended 
towards high F values. F determines the location of the 
peak in the matter power spectrum, and having extra red- 
shift information places tighter constraints on this - there 
is a degeneracy between F and the mean source redshift 
(z) such that a larger F would be compensated for by a 
smaller (z). Recall that F is not a fundamental quantity; 
in the limit of zero baryons, F = ^mh, which is modi- 
fied to F = a^/iexp (^-Ob (l + \/2/^/f]ni)), if f^b, the 
present baryon density, is accounted for in the transfer 
function (Sugiyama 1995). With our strong priors (hence 
fairly tight constraint on or a^), the error in F roughly 
translates into an error in /i, so redshift binning decreases 
the upper limit on h. 

In the ri,„ — (Tg plane we obtain the familiar "banana" 
shaped confidence regions, tightened with the inclusion of 
redshift binning. It is difficult to directly compare our con- 
straints to real surveys with different observational con- 
ditions; however, our confidence regions are roughly com- 
patible with those of van Waerbeke et al. (2001) allowing 
for these differences. 

To explore a wider range of parameter combinations, 
we employed the Fisher formalism since this allows one 
to easily obtain error estimates and investigate degenera- 
cies. We used the covariance matrices estimated from the 
first data set, again for 10 sub-fields as described above. 
Table [5] shows to what extent the errors on various pa- 
rameters are improved for — 2 compared with = 1. 
Note that these results depend on the cosmological model 
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and power spectrum of our fiducial model. Redshift bin- 
ning is particularly helpful when fewer strong priors are 
assumed, compared with the case when only one or two 
of the parameters are allowed to vary. In the case where 
rim, ^A, CTg and r are free, errors are a factor of roughly 
4 to 8 smaller when = 2. As we adopt more strong 
priors, redshift binning becomes progressively less benefi- 
cial. For example, if we consider parameter combinations 
where either fim or erg are assumed to be perfectly known, 
this breaks a strong degeneracy otherwise present; the de- 
crease in errors when Nz = 2 are therefore not so great as 
one might have anticipated. Another interesting trend is 
that the constraint of JIa seems to be most favourably af- 
fected by redshift binning, perhaps because it is important 
to the growth rate of structure at redshifts z ^ 1. 

How does Fisher analysis compare with the likelihood 
treatment? Fisher analysis should be seen as a way to esti- 
mate errors and investigate degeneracies, but does not pro- 
pose to reveal the detailed behaviour of confidence regions 
far from the fiducial model - for instance the asymmetric 
constraints in the F — 57,„ or —erg planes, or the "banana" 
shaped constraints in the ri„i — ug plane. Nevertheless, 
the differences in the constraints on the parameter pairs 
drawn from fim, dg and F are roughly consistent with the 
full likelihood treatment discussed above. 

Using the second data set consisting of 4 redshift bins, 
for a survey consisting of 1 sub-field, we obtained covari- 
ance matrices corresponding to various redshift cuts for 2 
and 3 redshift bins. Parameter constraints are again tight- 
ened going from iV^ = 1 to A^^ = 4. Considering again for 
example the case where all parameters are free, the most 
dramatic change is already seen going from Nz = 1 to 
A^2 = 2. In the case of A^z = 2, with the lowest redshift cut 
-^cut = 0.75 (2binsIII), parameter constraints are weakest. 
When higher redshift sources are isolated (zcut = 1-5, 2.25) 
the constraints are similar, and which is the better choice 
depends on the combination of parameters considered. For 
3 redshift bins, the combinations of cuts in redshift lead 
to very similar error estimates. 

At some point further sub-division into redshift bins 
does not lead to improved constraints on parameters (Hu 
1999). This limit must be determined for the survey and 
cosmological parameters in question. For simplicity con- 
sider the case of A^z = 2: and are correlated since 
sources in the upper redshift bin are also sensitive to lens- 
ing by structure at 2: < Zcut- A measure of this correlation 
is provided by a correlation parameter, p = ^lu/ (^lCu)*^^; 
also used by Hu (1999) in the context of power spec- 
tra, p = 1 and p = indicate complete correlation 
and lack of correlation respectively. For our simulations, 
0.83 < p < 0.93 for the range over which the correlation 
functions are calculated, with p = 0.87 taking the mean 
over all A9 bins. 

The improvement of the parameter constraints is 
clearly also a function of how the redshift intervals are 
set. This can already be seen in Table |31 where we calcu- 
lated different binnings with two and three redshift bins. 
We expect that there is an optimal way to bin the data. 



The intention of this paper, however, is to present a fast 
method for calculating the covariance matrix of the cosmic 
shear correlation estimator. This issue will be explored in 
a forthcoming publication. In practise, this question prob- 
ably does not arise anyway, because there one would take 
as many redshift bins as possible, their number being de- 
termined by the accuracy of the redshift estimator. 

As described above, one of our fiducial surveys is com- 
posed of 10 uncorrelated sub-fields. Another possibility 
is that a single contiguous patch of sky is targeted pri- 
marily for another science goal; such a survey might con- 
sist of 10 sub-fields drawn from the same large field. The 
sub-fields might in that case be selected so as to avoid 
bright stars, chip boundaries or defects. These sub-fields 
would be correlated to some extent, meaning that taking 
m rather than n(< m) sub-fields does not decrease the 
covariance by a factor of m/n, which would be the case if 
they were independent. The degree to which the sub-fields 
are correlated is accounted for directly in the covariance 
matrix. However, it is instructive to have an estimate of 
this, using the (ensemble average) covariance matrix for 
a survey composed of 10 correlated as opposed to 10 un- 
correlated sub-fields. The ratio of the diagonal elements of 
the correlated and uncorrelated geometry covariance ma- 
trices ranges between « 1.1 — 1.4, from the lower to upper 
redshift bin. Of course, the amount of correlation between 
sub-fields also depends on their geometry, so this has to 
be estimated for the survey in question. 

8. Conclusions 

In preparation for upcoming cosmic shear surveys, the 
main purpose of this paper was to demonstrate that it is 
possible to rapidly simulate surveys using a Monte Carlo 
method. This enables one to obtain accurate full covari- 
ance matrices for the shear two-point correlation function 
estimated from arbitrary survey geometries and with 
the sources binned in redshift. Averaging over many inde- 
pendent realisations enables us to take into account cosmic 
variance. As a first application, we estimated the extent 
to which redshift information for sources used in a cosmic 
shear analysis improves constraints on parameters derived 
from the estimated shear two-point correlation function 

A likelihood analysis in the fim — cg, f2„i — T and F — ag 
planes shows that separating the sources into two redshift 
bins enables tighter constraints to be placed on parame- 
ters. Considering a wider range of parameter combinations 
in the context of a Fisher analysis reveals that redshift in- 
formation is particularly advantageous in cases where few 
strong priors are assumed. When ri,„, fS: T and rt\ are 
free parameters, having 2 (4) redshift bins tightens errors 
on parameters by a factor of ~ 4 — 8 (~ 5 — 10). Most im- 
provement on error estimates occurs going from Nz = 1 
to A^z = 2. In general, for the combinations of free and 
fixed parameters explored, Oa seems to benefit most from 
redshift binning. 



Simon, King & Schneider: 3D correlations of cosmic shear 



13 



One might ask why cosmic shear is of interest, in the 
hght of the recent WMAP results (e.g. Bennett et al. 
2003), which suggest that the era of precision cosmology 
is already upon us; there are several facets to consider. 
Cosmic shear has the power to break degeneracies inher- 
ent to CMB data (e.g. Hu & Tegmark 1999), for instance 
the angular diameter distance degeneracy (e.g. Efstathiou 
& Bond 1999). It also provides a completely independent 
cross-check of cosmological parameters, based on equally 
well understood but different physical principles. Besides 
this, as is the case with the CMB, the interpretation of 
results requires no assumption about the bias between lu- 
minous tracers and the underlying dark matter distribu- 
tion which plagues, for example, galaxy redshift surveys. 
In addition, with redshift estimates for sources in upcom- 
ing large cosmic shear surveys, lensing has the potential 
to see beyond the radially projected convergence power 
spectrum, becoming sensitive to the evolution and growth 
of structure in the Universe. 
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Appendix A: Switching from a finer to a coarser 
redshift binning 

Here we show how the auto- and cross-correlations of the 
cosmic shear from a finer redshift binning are related to 
the auto- and cross-correlations obtained from a coarser 
redshift binning (by combining the finer bins). The re- 
lations of this section ensure that we only have to make 
simulated data of the finest redshift binning, since the cor- 
responding data with less redshift information can always 
be related to the SD-correlations of the cosmic shear of 
this case. 

In a first step, we turn to the auto-correlation ^± of 
a new redshift bin, neglecting for a moment the cross- 
correlations to the other new redshift bins. ^± is according 
to the Eqs. ^ and (0 a function that linearly depends on 

(w) 



^± (0) 



4c4 

du jw) 

— dw Jo.4 [£0) Ps 

27r Jq (w) 



(A.l) 



/H 



If we split the redshift distribution p (w) of the source 
galaxies into disjunct parts, like 



g« (w) 



p{w) for w e [wi-i,Wi 
else 



with Wi{zi) in the sense of Fig.^ we expand with the no- 
tation of Eqs. Q the function W in the following manner 



W"^ (w) = ^ iMUjW'^'^ (w) {w) 



(A.3) 



As the g*^*-* defined here are not normalised, but the cor- 
responding redshift distributions p^*^ in the definition for 
M^(') (w) are, we introduce the normalisation factors 



dw g'^'-' (w) 



(A.4) 



telling us what fraction of the distribution inside the 
new bin is contained in its subdivisions q^'^\ Translating 
Eq. IjA.Sp to the power spectra gives, leaving out the ar- 
guments in £ 



(A.5) 



i>j 



Similarly, we get for the cosmic shear auto-correlation 
^ (actually for all linear functions of Pk''''; hence the 
dropped index "±"): 



■J S,4J : 



(A.6) 



where S^a are the auto-correlations for the sub-bins, S^ij 
are the cross-correlations of the cosmic shear between the 
sub-bins. This equation tells us, therefore, how we have 
to combine the cosmic shear correlations of the sub-bins 
when we are switching from a finer to a coarser redshift 
binning of the data. 

What about the cross-correlations between the new 
redshift bins if we decide to switch to a binning with more 
than one redshift bin? This case is treated like the fore- 
going one, except that it is slightly more general. Assume 
we focus on two new redshift bins k and / consisting of 
data from a finer redshift binning: 

M = E 'i' (^) (w) = E (^) ' (A-7) 



where pk is the redshift distribution inside the new bin k 

(i) 

and pi the same for the new bin I. Bin k combines q^, and 

the bin I combines from a finer binning, respectively. 
Using the same arguments as before, we obtain the follow- 
ing relation between the cosmic shear cross-correlation be- 
tween the new redshift bins £,ki , and the cross-correlations 
4ki) 



between their components: 



iki - E 



(fc) {i)Aki) 



(A.8) 



n,-'^-' and n[^^ are the normalisations for the sub-bins. This 

(A 2) ■' n I, 

^ ■ ' equation is, of course, the generalisation of Eq. ljA.6|) . 
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